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ABSTRACT 


This  is  a  continuation  of  a  heat-transfer  investigation 
performed  as  In-House  Laboratory  Independent  Research  to  de¬ 
termine  the  effect  of  interface  thermal -res i stance  on  multi¬ 
layer  gun  barrel  radial  temperature  distributions.  In  pre¬ 
vious  Technical  Reports  69-121  and  70-155,  published  under 
the  same  title,  the  task  of  establishing  the  feasibility  of 
induced  interface-resistance  was  discussed. 

During  this  reporting  period,  heat-transfer  analysis  and 
experimental  correlation  were  continued.  An  experimental 
M60  gun  barrel  was  fabricated  and  tested.  The  reduced  outside 
barrel  temperatures,  of  this  firing  test,  demonstrated  the 
effect  of  interface  thermal  resistance.  Also,  a  significant 
effort  was  directed  toward  establishing  fabrication  techniques 
for  a  full-length  thermal  interface  barrel.  As  a  result  of 
these  efforts,  a  patent  application  has  been  submitted  by 
the  Research  Directorate,  Weapons  Laboratory  at  Rock  Island 
to  cover  this  concept. 
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OBJECTIVE 


The  objective  of  this  analysis  was  to  analytically  and 
experimentally  determine  the  effect  of  interface  thermal 
resistance  on  multilayer  gun  barrel  temperature  distri¬ 
butions.  Specifically,  analysis  was  to  be  conducted  to 
establish  the  design  criteria  for  an  experimental  test 
barrel.  Then,  based  on  experimental  correlation  and  ad¬ 
ditional  analysis,  design  requirements  were  to  be  deter¬ 
mined  for  a  full-length  modified  gun  barrel  with  consider¬ 
ation  of  the  interfacial  material  requirements  and  forming 
processes. 

INTRODUCTION 


Consistent  with  an  overall  effort  to  increase  weapon- 
firepower,  this  heat  transfer  analysis  performed  by  the 
Research  Directorate,  Weapons  Laboratory  at  Rock  Island  has 
been  directed  toward  reducing  gun  barrel  heat-transfer. 
Previous  investigations  under  this  subject,  have  been  di¬ 
rected  toward  the  preliminary  task  of  analytically  and  ex¬ 
perimentally  substantiating  the  Influence  of  Induced-inter- 
face  thermal  resistance,  and  of  measuring  Interface  thermal 
resistance  for  various  Interfacial  conditions. 

The  current  Investigation  was  coupled  with  previous 
analysis  to  establish  design  criteria  for  a  full-length 
composite  gun  barrel  with  Induced  interface  resistance.  An 
N60,  7.62mm  gun  barrel  was  designed,  fabricated  and  both 
analytically  and  experimentally  evaluated. 

A  comprehensive  study  was  undertaken  to  evaluate  various 
fabrication  techniques  for  gun  barrels  with  full-length 
thermal-resistance  Interfaces.  As  part  of  this  study,  a  ma¬ 
terial  survey  was  performed  to  establish  the  most  promising 
Interface  materials.  Parameters  considered  In  this  survey 
included  thermal  conductivity,  material  stress  properties, 
and  thermal  coefficient  of  expansion.  The  various  barrel¬ 
forming  techniques  evaluated  Included  co-extrusion,  swaging, 
and  shrink-fit.  Guidelines  were  established  for  future 
fabrication  of  full-length  modified  gun  barrels. 

Acknowledgment  Is  extended  to  the  Metallic  Materials 
Application  Team  and  to  the  Process  Technology  Application 
Team,  both  of  the  Research  Directorate,  Weapons  Laboratory 
at  Rock  Island  for  their  contributions  and  guidance  in  the 
establishment  of  multilayer  gun-barrel-forming  criteria  and 
In  the  application  of  interface  material  coatings. 
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ANALYTICAL  AND  EXPERIMENTAL  INVESTIGAT 1 0 N 


Heat  Transfer  Analysis 

As  part  of  the  effort  to  establish  design  criteria  for 
a  full-length  thermal  Interface  gun  barrel,  analysis  was 
performed  on  an  M60  gun  barrel  to  calculate  temperature 
distribution  as  a  function  of  Induced  Interface  thermal  re¬ 
sistance.  A  gun  barrel  with  a  configuration,  as  shown  In 
Figure  1,  was  analyzed  to  predict  radial  temperature  dis¬ 
tributions  In  the  modified  section  of  the  gun  barrel. 

Experimental  time-temperature  history  data  for  a  stan¬ 
dard  M60  gun  barrel  was  used  to  calculate  effective  pro¬ 
pellant  gas  convection  coefficients  and  gas  temperatures. 
For  the  condition  In  which  all  heat  Is  dissipated  In  the 
gun  barrel  or  In  which  heat  Input  by  convection  Is  equal  to 
the  change  In  barrel  Internal  energy,  the  governing  expres¬ 
sion  Is: 

hgAtYT]  *  pVrAi»  35 


where  hg  *  propellant  gas  convection  coefficient 


Tg  "  propellant  gas  temperature 


A  *  bore  surface  area 


Am  *  mean  surface  area 
m 


T  »  gun  barrel  temperature 


Ar  *  wall  thickness 


p  *  gun  barrel  material  density 


specific  heat  of  gun  barrel  material 


rate  of  change  of  gun  barrel  temperatures 
with  respect  to  time. 
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FIGURE  1 


As  the  outside  barrel  temperature  increases,  dissipation  to 
the  surrounding  environment  must  be  accounted  for  in  the 
equation.  Heat  transfer  by  radiation  and  convection  is  ad¬ 
ded  to  the  equation  and  the  solution  for  h  and  T  follows 
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by  the  same  procedure.  A  typical  computer  program,  based 
on  the  equation  given  above,  for  the  calculation  of  effec¬ 
tive  propellant  gas  convection  coefficients  and  temperatures 
is  given  in  Appendix  A, 

Knowing  effective  propellant  gas  temperatures  and  con¬ 
vection  coefficients,  radial  gun  barrel  temperatures  (as  a 
function  of  rounds  fired)  can  be  calculated  for  composite 
gun  barrels  with  i.'.luced  interface  resistance.  A  computer 
program  in  which  a  Crank-Ni col  son  algorithm  is  used  was 
employed  for  this  purpose.  A  copy  of  this  program  is  given 
in  Appendix  B,  Results  of  these  calculations  are  shown  in 
Figure  4. 

Based  on  the  induced  interface  thermal  resistance  re¬ 
quirements  determined  by  this  analysis,  an  M6Q  gun  barrel 
was  modified  by  the  insertion  of  a  Zirconium  Oxide  (Zr02) 
coated  section  as  shown  in  Figure  1.  The  various  fabri¬ 
cation  and  assembly  drawings  for  this  modified  gun  barrel 
are  shown  in  Appendix  C.  Detailed  steps  involved  in  the 
fabrication  of  the  insert  section  and  in  the  application 
of  the  Zr02  coating  are  shown  in  Figure  2,  This  was  the 
first  attempt  to  shrink-fit  an  outer  sleeve  over  a  Zr02 
coating.  Results  of  this  effort  were  applied  in  the  fabri¬ 
cation  investigation  of  a  full-length  modified  gun  barrel 
discussed  later  in  this  report 

Experimental  Results 

The  modified  M60  gun  barrel  of  Figure  1,  along  with  a 
standard  barrel,  was  instrumented  with  thermocouples  for  a 
live-firing  experiment  to  compare  barrel  temperatures  and 
as  an  experimental  correlation  of  the  theoretical  analysis. 
This  experiment  consisted  of  continuous  automatic  fire  of 
400  rounds  at  a  firing  rate  of  600  rounds  per  minute. 

Results  of  this  firing  experiment  are  presented  in 
Figures  3  and  4.  External  gun  barrel  temperatures  are  plotted 
for  two  axial  locations,  11.2  and  12,2  inches  from  the  breech, 
both  located  within  the  modified  section  of  the  gun  barrel. 
Both  standard  and  modified  measured  gun  barrel  temperatures 
at  the  12.2-inch  location  are  given  in  Figure  3;  whereas, 
both  standard  and  modified  measured  gun  barrel  temperatures 
at  the  11.2  location,  along  with  calculated  temperatures 
for  this  location  are  compared  in  Figure  4, 
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Zirconium  Oxide  Applied  Over  Coated  Liner  Ground  to  Fit  C140299 

Bond  Coat  OD  =  .588" 

OD  =  .612" 


FIGURE  2 


1200 


ROUNDS  FIRED 


MEASURED  M60  EXTERNAL  TEMPERATURE 


ROUNDS  FIRED 


The  modified  section  operated  at  approximately  100°F 
cooler  after  400  rounds  than  the  same  location  on  a  stan¬ 
dard  gun  barrel.  Also,  calculated  temperatures,  shown  in 
Figure  4,  are  lower  than  those  measured.  This  can  be 
accounted  for  because  the  axial  heat  conduction  into  the 
specimen  was  neglected  in  the  calculations.  A  signifi¬ 
cant  amount  of  heat  is  transferred  from  the  higher  tempera¬ 
ture,  unmodified  section  into  the  modified  section.  Hence, 
a  full-length  modified  gun  barrel  would  operate  at  lower 
temperatures  than  measured,  that  is,  at  temperatures  closer 
to  those  of  the  analytical  calculations. 

INVESTIGATION  OF  FABRICATION  TECHNIQUES  FOR  MULTILAYER  GUN 
BARRELS  WITH  INDUCED  INTERFACE  THERMAL  RESISTANCE 


The  investigation  of  fabrication  techniques  for  multi¬ 
layer  gun  barrels  with  induced  thermal  resistance  was  sub¬ 
divided  into  two  interrelated  subtasks  -  material  selection 
and  forming  processes.  Material  properties  will  be  par¬ 
tially  altered  by  the  forming  process;  furthermore,  because 
of  the  limited  number  of  low  thermal  conductivity  materials, 
the  number  of  forming  processes  available  is  limited  by 
material  properties. 

Interface  Material  Selection 


Various  ceramics  and  other  materials  have  been  investi¬ 
gated  for  application  in  the  fabrication  of  a  full-length 
thermal  interface  gun  barrel.  Hi gh-temperature  resistant 
materials  with  the  low-thermal  conductivity  (less  than 
1  BTU/hr  ft  °F)  required  for  a  thermal  barrier  typically 
exhibit  low  tensile  and  shock  strength  properties.  Nor¬ 
mally,  plastic  deformation  is  absent  prior  to  mechanical 
failure.  With  these  property  limitations  in  mind,  two  de¬ 
sign  approaches  have  been  taken  to  obtain  a  mechanically 
stable  barrier.  First,  the  material  can  be  a  load-carrying 
member,  solid  configuration,  capable  of  withstanding  hoop 
stresses,  thermal  stresses,  and  axial  shear  stresses  present 
during  normal  firing.  Secondly,  the  material  can  be  required 
to  transmit  loads  to  the  outer  sleeve  that  functions  as  a 
support  for  the  entire  assembly. 

Based  on  these  design  requirements,  thermal -materi al 
selection  criteria  have  been  established  as  follows: 

1.  Material  must  have  low-thermal  conductivity  (less 
than  1  BTU/hr  ft  °F). 
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2.  Material  must  be  chemically  and  metal  1 urgical ly 
stable,  up  to  3000°F.  That  is,  no  chemical  or  metallurgical 
reactions  occur,  such  as  material  phase  change  or  carbiding 
of  other  composite  materials. 

3.  Material  should  have  high-thermal  shock  resistance 
to  preclude  material  breakdown  during  repeated  thermal 
cycling  (primarily  for  the  solid  configuration  assemblies). 

4.  Interfacial  material  must  have  a  linear  coefficient 
of  thermal  expansion  (CTE)  similar  to  that  of  the  barrel 
materials  (Cr-Mo-V  steel,  a  typical  barrel  material  having 

a  CTE  of  7.8  X  1Q"S  in/in  °F).  Any  CTE  mismatch  can  result 
in  failure  in  both  the  axial  and  the  radial  directions  be¬ 
cause  of  assembly  separation  during  thermal  cycling. 

5.  Interface  material  must  be  capable  of  sustaining 
or  transmitting  all  modes  of  barrel  loading. 

The  two  most  promising  interface  materials  that  closely 
satisfy  these  criteria  are  zirconium  oxide  (ZrO,,)  and  aniso¬ 
tropic  boron  pyrolytic  graphite  (bPG)„  Because  of  the  dif¬ 
ferent  nature  of  these  two  materials,  they  will  be  discussed 
separately. 

Zirconium  Oxide 


Zirconium  oxide  (Zr02),  dependent  upon  the  manufacture 
of  the  basic  oxide  and  the  coating  process,  varies  in  me¬ 
chanical  and  thermal  properties.  Material  density  (percent¬ 
age  of  porosity),  degree  of  stabilizing  agent  employed 
(MgO,  CaO  prevent  material  phase-change  throughout  the 
operational  temperature  range),  and  degree  of  metallic 
alloying  are  the  primary  causes  of  these  property  variances. 
Zr02  may  be  found  in  various  commercial  forms,  some  of 
which  are  as  follows: 

1.  Powders  -  applicable  for  flame-spraying. 

2.  Free-standing  ceramic  bodies  (for  example,  extruded 
or  pressed  parts). 

3.  Tapes,  papers,  and  fabrics. 

A  tabulation  of  typical  thermal  and  mechanical  proper¬ 
ties  for  a  sintered  ZrO?  (flame  spray,  partial  stabilization) 
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is  as  foil ows : 


1.  Thermal  conductivity  -  541  BTU/hr  ft  °F 

2.  Oensity  -  374  L bm / F t 3 

3.  Linear  coefficient  of  thermal  expansion 
5.6  X  10"6  in/in  °F  (RT  to  2200°F) 

4 .  Specific  heat  0.13-0.18  BTU/Lb  °F  at  200  to  1950°F 

m 

5.  Ultimate  tensile  strength  (psi) 

21,120  at  RT 
1  5,000  1  6  2  5 0  F 

13,230  1 886°F 

12  ,000  21  92°  F 

6.  Ultimate  compressive  strength  (psi) 

303,000  at  RT 
100,000  2370°  F 

For  a  load-sustaining  assembly  (solid  material  config¬ 
uration),  the  ceramic  material  must  maintain  its  structural 
integrity.  Therefore,  the  material  mechanical  properties 
must  be  sufficient  to  withstand  all  stress  modes.  Initial 
stress  calculations  with  a  homogeneous  barrel  assembly  in¬ 
dicate  that  material  failure  will  occur  unless  a  precom- 
pressive  stressing  of  the  interface  material  can  be  obtained. 
For  example,  a  gas  pressure  of  65,000  psi  will  cause  a  tensile 
stress  at  the  interface  in  excess  of  25,000  psi  that  exceeds 
the  ultimate  tensile  strength  of  ZrOt  Also,  uneven  axial 
loading  can  be  a  source  of  interface-material  cracking.  Ex¬ 
perimental  verification  of  assembly  structural  integrity 
must  be  performed  to  determine  the  significance  of  cracking.  • 

In  a  load  transmission  assembly,  l' 0,  is  required  during 
forming  to  compress  to  a  highly  compacted,  dense  powder. 

During  actual  firing,  no  further  compaction  or  shifting  of 
the  material  is  allowed  The  assembly  procedure  necessary 
to  obtain  maximum  compactions  and  proper  natural  flow  will 
be  determined  experimentally  in  the  future. 

Anisotropic  Boron  Pyrolytic  Graphite  (bPS) 

Boron  pyrolytic  graphite  (bPG)  exhibits  strong  aniso¬ 
tropy  (deferent  property  values  along  different  axes)  In 
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thermal  and  mechanical  properties.  This  material  is  manu¬ 
factured  by  a  high-temperature,  vapor-deposition  process; 
therefore,  it  occurs  only  as  a  free-standing  body  con¬ 
figuration.  The  diagram  below  indicates  the  primary  prop¬ 
erty  directions  for  a  radial  element; 


c 


b 


Radial  Element 


Anisotropy  Is  exhibited  In  only  two  directions  of  a  cylin¬ 
drical  coordinate  system.  The  "c"  property  direction  lies 
on  the  radial  component  and  the  "ab!l  property  direction  Is 
a  surface  comprising  the  circumference  and  the  axial  axis. 
For  use  as  an  Interface  material,  the  "cM  direction  Is  the 
controlling  direction.  The  alloying  of  basic  pyrolytic 
graphite  with  boron  (0  to  1.2%)  results  In  alteration  of 
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material  properties  in  both  the  "ab"  and  "c"  directions. 

A  percentage  increase  in  boron  causes  a  decrease  in  thermal 
conductivity,  a  major  advantage  for  this  application. 

A  material  property  compilation  is  presented  in  Table  I 
for  a  zero  percentage  boron  content: 

TABLE  I 


Direction 


Property 

£ 

ab 

Units 

Ultimata  Tensile 

1  ,500 

14,000 

psi 

Strength  ( AT ) 

Ultimate  Flexural 

1  ,500 

21 ,000 

psl 

Strength  (RT) 

Ultimate  Compressive 

68,000 

14,000 

psi 

Strength  (RT) 

Thermal  Conductivity  (RT) 

1 .00 

200 

BTU/hr  °F 

(2500°F) 

0.66 

58 

Specific  Heat  (RT) 

0.26 

BTU/Lbn,  ®F 

Linear  Coefficient 

10.8 

0.2 

In/In  ®F  X 

of  Thermal  Expansion 


Thermal  cond'«ct1  vi ty  data  are  given  in  Table  II  for  the  "c" 
direction  as  a  function  of  the  percentage  of  boron  alloying: 


TABLE  II 


Per  Cent,  Boron 
1.2 
075 
0.35 
0 


Thermal  Conductivity 
BTU/hr  ft  ®F 

v  0  36 

045 

0  54 

1  00 
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Mechanical  properties  are  unavailable  for  boron  alloying 
Theory  indicates  that  alloying  will  increase  strength 
properties. 

Application  of  bPG  as  an  interface  materia!  is  feasible 
for  a  solid  configuration.  Some  work  has  been  done  with 
co-extrusion  of  isotropic  graphite  materials.  The  diffi- 
culty  encountered  was  that  isotropic  graphite  did  not  de¬ 
form  plastically  during  extrusion;  therefore,  composite 
failure  occurred. 


A  deterrent  factor  in  applying  bPG  as  an  interface 
material  Involves  the  mismatch  of  linear  coeffic’ents  of 
thermal  expansion  (CTE).  A  comparison  of  bPG  jaterial  with 
typical  barrel  materials  indicates  that,  both  in  the  c 
direction  and  on  "ab"  surface,  the  CTE  vanes  significantly 
enough  to  cause  composite  failure  because  of  material  sepa¬ 
ration  during  thermal  cycling  Also,  b*ndmg  during  shrink- 
fit  operations  could  occur  due  to  expansion  in  the  c  di¬ 
rection.  Further  analysis  and  experimentation  will  have 
to  be  undertaken  to  confirm  and  'eso*ve  these  problems 


during 

carbon 


The  problem  of  carbonizat ’on  of  ba-rel  steel 
thermal  cycling  was  investigated  Nuc’eation  of 
atoms  along  the  steel  gram  boundaries  could  initiate 
microscopic  cracking  due  to  mater»a’  hardening,  but  the 
feasibility  of  composite  failure  due  to  carburizing 
termlned  minimal 


was  de- 


Forming  Processes 

The  primary  consideration  m  the  forming  process  is  that 
of  the  structural  integrity  of  the  composite  gun  tube. 

Basic  design  considerations  used  as  guidelines  for  this  in¬ 
vestigation  are  listed  below? 

1.  Liner  and  sleeve  assemblies  must  be  the  primary 
load-carrying  composite  members 

2.  Feasibility  must  be  determined  fo»  application  of 
precompress ive  loading  at  the  interface  to  offset  tensile 
stresses  generated  by  pressure  and  thermal  stressing 

3.  Significant  design  changes  to  the  overall  weapons 
system  must  be  kept  to  a  minimum 

4.  Bore  liner  material  must  withstand  the  higher  bore 
temperatures  engendered  by  application  of  the  therma 
barrier. 
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Three  principal  fabrication  processes  have  been  In¬ 
vestigated.  These  are  as  follows: 

1.  Shrink-fit  assembly  (one-  or  two-piece  approach). 

2.  Swaging  of  outer  sleeve  over  inner  liner  and 
interface  material . 

3.  Hot  coextrusion  of  entire  three-piece  assembly. 

All  these  processes  require  three  major  components,  an 
inner  liner,  a  thermal  Interface  resistance  material,  and 
an  outer  retaining  sleeve.  Each  of  these  processes  will 
be  discussed  In  the  following  sections. 

Shrink  Fit 

Current  manufacturing  procedures  limit  the  shrinkage 
length  to  approximately  10  Inches.  By  an  Increase  In  the 
assembly  rate,  an  Increase  In  the  temperature  difference 
between  sections,  and  an  Increase  In  the  clearance  toler¬ 
ances  (a  reduction  In  final  Interference  between  sections), 
longer  length  fits  can  be  achieved.  By  controlled  procedures, 
sections  have  been  assembled  In  excess  of  30  Inches.  Initial 
design  work  on  an  M60  gun  barrel  Indicates  that  a  one-piece 
shrink-fit  can  be  accomplished. 

The  primary  advantages  of  a  shrink-fit  assembly  are 
listed  below: 

1.  Solid  interface  sections  may  be  employed,  either  by 
flame  spraying  the  Interface  material  on  the  liner,  by  use 
of  free-standing  ceramic  cylinders,  or  by  use  of  ceramic 
tapes  or  fabrics. 

2  A  precompresslve  load  may  be  applied  (dependent 
upon  the  amount  of  Interference  desired)  to  the  ceramic  In¬ 
terface,  thus  the  initial  stress  characteristics  are  enhanced. 

With  a  solid  interface  section  configuration,  the  bond¬ 
ing  mechanism  between  liner  and  Interface  material,  and  like¬ 
wise  between  Interface  material  and  outer  sleeve,  should  be 
considered  Bonding  may  be  established  through  mechanical 
means  such  as  roughened  surfaces,  cements,  or  flame-spray 
coats. 

Initial  calculations  performed  for  an  N60  gun  barrel 
show  that  a  precompresslve  load  of  59,000  psl  can  be  devel¬ 
oped  for  a  1,5-mil  Interference.  An  expansion  of  2.5-m11$ 
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can  be  achieved  with  a  temperature  difference  of  1000°F. 
Hence,  by  an  Increase  In  assembly  rates,  a  longer  shrink- 
fit  apparently  can  be  achieved. 

Swaging 

Swaging  has  a  direct  application  for  the  forming  of  a 
composite  material  barrel  blank.  The  process  Is  relatively 
Inexpensive  when  compared  with  other  processes.  .Major  prob¬ 
lems  are  generated  In  the  swaging  of  brittle  materials. 

Rotary  swaging  will  cause,  as  a  function  of  radial  angle  * 
and  axial  speed  of  the  workpiece,  an  uneven  reduction  In 
material  cross-section.  The  cyclical  nature  (a  fluctuating 
angular  radial  force)  of  the  swaging  process  and  the  low 
Impact-resistance  of  ceramic  materials  will  cause  pulveri¬ 
zation  of  the  ceramic  Interface  during  assembly*  thi s  results 
In  the  ejection  of  the  ceramic  material  from  the  assembly. 
Also,  Interface  material  will  be  unevenly  distributed  around 
the  circumference  of  the  assembly:  Even  though  solutions  to 
the  many  problems  associated  with' the  swaging  process  have 
not  yet  been  determined,  this  process  is. believed  to  have 
excellent  potential  In  view  of  the  low  process  cost  Involved. 

Hot  Extrusion 

Under  normal  practice,  when  two  or  more  materials  are 
extruded  together,  they  become  metal lurgically  bonded  at 
the  interface.  In  the  case  of  ceramics  and  metals,  this 
action  does  not  occur.  Ap  end  configuration  comprises  an 
Inner  liner,  a  densely  compacted  oxide  powder,  and  an  outer 
sleeve.  Because  of  toe  differences  in  hot  working  tempera¬ 
tures  (Cr-Mo-V  steel  frod  19pO°F  to  2100’F  and  ZrO*  approx¬ 
imately  4800°F) ,  plastic  deformation  o'  the  ceramic  material 
Is  Impossible  without  exceeding  Its  ultimate  shear  strength. 
The  resulting  Interface  transmits  hoop  stresses  and  thermal 
stresses  to  the  strength-bearing  members.  However,  axial 
shear-stresses  due  to  projectile  frictional  forces  can 
cause  assembly  failure 

6 

A  number  of  approaches  fo-  the  solution  to  the  problem 
of  axial  shear  failures  are  as  fQllows:( 

1  Electron  beam-welding  of  retaking  rings  over  the 
muzzle  and  breech  ends  of  the  gun  barrel  This  approach  Is 
applicable  to  all  processes  under  discussion, 

2.  Another  approach,  applicable  to  hot  extrusion,  would 
Involve  the  use  of  liner  supports  These  supports  would  be 
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attached  periodically  along  the  length  of  the  liner  (similar 
to  cooling  fins).  Ceramic  material  would  be  used  to  make 
up  the  remainder  of  the  annulus.  The  liner,  liner  support, 
ceramic,  and  sleeve  section  would  then  be  extruded.  The 
resulting  configuration  would  consist  of  a  liner  support, 
metal! urgical ly  bonded  to  both  the  liner  and  the  sleeve 
assembly,  with  the  interface  material  distributed  between 
the  liner  support  as  a  high  density  powder. 

3.  The  interface  material  can  be  alloyed  with  free 
metallic  atoms,  copper  for  example.  The  percentage  of 
alloying  can  be  varied  from  10  to  20  per  cent  of  the  inter¬ 
face  material.  During  extrusion,  the  alloyed  material  will 
flow  and  form  metallurgical  bonds  with  the  liner  and  the 
sleeve  assemblies.  The  resultant  configuration  would  con¬ 
sist  of  a  solid  interface  material  that  could  withstand 
axial  shear. 

FUTURE  EFFORTS 


Presently,  two  gun  barrels  are  being  evaluated  for  fab¬ 
rication  of  a  full-length  thermal  resistance  interface.  The 
two  gun  barrels  are  an  M60  7.62mm  and  a  MK1 1  Mod.  5, 20mm. 
Preliminary  design  of  an  M60  gun  tube  with  a  shrink-fit 
assembly  (30  mil  coating  of  Zr02,  high-density  flame  spray) 
is  nearing  completion.  A  major  redesign  of  the  breech  sec¬ 
tion  is  necessarily  unavoidable  to  inclose  a  full-length 
interface  barrel.  An  extensive  test  program  will  be  initi¬ 
ated  to  measure  temperature  profiles  and  erosion  wear  as  a 
consequence  of  high  bore-temperatures . 

A  coextrusion  process  is  under  investigation  for  placing 
a  30-  to  35-mil  Zr0?  interface  into  a  MK1 1  Mod.  5  gun  barrel. 
This  investigation  includes  selection  of  the  liner  material 
and  the  outer  sleeve  material  consistent  with  the  coextru¬ 
sion-process  and  the  operational  temperatures.  Three  com¬ 
plete  MK1 1  gun  barrels  are  contemplated  for  fabrication  and 
testing. 

SUMMARY  AND  CONCLUSIONS 


A  second  live-firing  experiment  was  completed  using  an 
M60  gun  barrel  with  a  Zr02  thermal  interface  as  the  test 
vehicle.  The  test  section  for  this  experiment  was  located 
near  the  muzzle  end.  Good  correlation  between  theory  and 
experimentation  was  obtained.  A  significant  number  of 
rounds  were  fired  with  no  apparent  deterioration  of  the 
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modified  section.  Experimental  results  further  substantiate 
the  temperature  reduction  for  the  outer  sleeve  of  the  modi¬ 
fied  section.  The  significance  of  this  temperature  reduction 
becomes  more  apparent  when  axial  conduction  into  the  test 
section  from  the  unmodified  portions  of  the  gun  barrel  are 
taken  into  account. 

Design  criteria  have  been  established  for  two  full-length 
thermal  interface  gun  barrels.  The  7.62mm,  M60  gun  barrel 
will  consist  of  a  shrink-fitted  outer  sleeve  over  a  coated 
liner.  The  20mm  design,  however,  will  be  based  on  a  coex¬ 
truded  composite  gun  barrel  with  ZrO?  positioned  between  a 
high-temperature  liner  material  and  an  outer  sleeve  of 
Cr-Mo-V  steel.  Future  efforts  by  the  Research  Directorate, 
Weapons  Laboratory  at  Rock  Island  will  involve  the  fabri¬ 
cation  and  testing  of  these  full-length  modified  gun  barrels. 
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//PRINT  JOB  (0174093, 1>, 'GEO  STILES'  "  '  JOB  "221 


C  PROGRAM  TO  DETERMINE  VALUES  POR  HG  £  TG  2/1B/71 

DIMENSION  9B(20)j  CC  (  420)^  DO  (420)  ,  _  _ HBAR  (  S  ) ,  TW  (  B  ) ,  DTWDT  1  r> ) _ 

l»  TB(5) ,  T0<5), T"n  f  20 ) , ~  CPl(20»V"  yT(30I 
COMMON  /BLK2/  X(30),  Yl(30),  Y2I30),  DTDT( 30) *  YYI30) 

TW  -  AVERAGE  WALL  TEMPERATURE 
RHO  -  OENSITY 
CP  -  SPECIFIC  HEAT 

HBAR  -  EFFECTIVE  CONVECTION  COEFFICIENT 
HBAR  =  XK  1  *  TW**2  +XK2  *  TW  *XK3 
DR  -  DELTA  R 

T A  -  AMBIENT  TEMPERATURE 
THETA  -  TIME 

R  I  IS  BORE  R  ADI  US  _  _  _ _ 

.  RO  IS"  OUTSIDE  RADIUS  .  . . 

RATIO  IS  FROM  1  THRU  KU/Ri 

IKKMAX  REPRESENTS  NUMBER  OF  -TIME  VS  TEMP-  DATA  CHANGES 
IKKMAX  =  I 

DO  111  IKK  =  1, IKKMAX 
C  N  REPRESENTS  NUMBER  OF  ORDERED  PAIRS. 

REAO  1,  N,  (X(H,  Yl(f)-,"i  =1,N) 

1  FORMAT! I5/I2F10.0) ) 

C  M  IS  THE  NUMBER  OF  CONSTANTS  NEEDED  FOR  DEGREE  OF  POLYN.  DESIRED. 

M  =  3 

C  NC  IS  THE  NUMBER  OF  PASSES  THROUGH  LEAST  SUUARES  SUBROUTINE,  AT  THIS 

C  POINT  IS  USED  ONLY  TO  MATCH  TI MF  VS.  TFMPERATURE. 

"""  "n‘c'“»T . * . . . . .  .  ~~ . ~ . ‘ 

CALL  LSTSQI M»N»NC ) 

DO  3  IJ  *  l,  N 
3  YI( IJ)  =  YY ( IJ) 

READ  5, IRATIOfKI ,  RO,  RHO,  CP,  XK1,  XK2,  DR,  TA,  XK3 
5  FORMAT ( I5/( 8F10.0) ) 

DR  =  RO  -  RI 


C  IRATIO  IS  THE  NUMBER  OF  RATIO  TESTS  PER  EACH  IKK 


DO  110  I J I  =  1,  IRATIO 

READ  8,  RATIO 

8  FORMAT ( F 10. 0 ) 

RO  =  RATIO  *  RI 

C 

11  =  0 

C 

N  IS  THE  NUMBER 
NN  =  N  -  1 

OF  DATA  SETS. 

c 

NN  IN  NEXT  CA  <0 
DO  100  IJ  =  l.NN 

MUST  BE  CORRECT  NUMBER  OF  DATA  SETS  READ  IN  -  ONE  *  * 

ram  =  ykiji 

DTWDT  ( 1 )  =  DTDT  (  IJ)  *  3  60_U,j) 


YUU  +  l) 


o  o  o 


i  DTWDTI2)  =  DTDTIIJ  ♦  1)  *  3600.0 

C  __  _ _ _ _ 

PRINT  20,  RATIO,  TOI 1) , TO ( 2 ) , D  T  WO  Till ,  DTWDT I  2 ) , RHO,  CP.XKl, 

_ 1  XK2.  DR.  TA,  XK3 ,  Rl,  RO. _ 

20  F0RMAT(10H  RO  /  RI  =,F12.4,  10X,  8H  TOIi)  =,F12.4, 

1  10Xj_  8H  TO ( 2 )  = , F 12 . 4/7H  Pr(l)=,FI2.4.  IPX,  6HDT(2>=,F12.4,  IPX, 

2  6H  RHO  =,FL2.4,  IOX,  5H  CP  =.F12.4/  5H  KI  =,E12.4,5X,5H  K2  =, 

_  3  E 12.4,  5X,  5H  DR  =,  E12.4,5X,5H  TA  =,F8.l/I5X,  6H  XK3  =,  EI2.&, 

4  10X.14H  BORE  RADIUS  =,EI2.6,  IOX,  17H  OUTSIDE  RADIUS  =,  E12.6//1 

_ 00  30  I  =  I,  2 _ 

HBAR(I)  =  XK1  *  TO(I)  **  2  *  XK2  *  TOII)  +  XK3 

.  PRINT  25,1, HBAR (  1  ) _ _ 

25  FORMAT!  8H  HBAR (  ,I2,5H  )  =  ,E12.5) 

30  CONTINUE  _  _ 

AB  =  2.  *  RI 

_ AO  =  2.  *  RO _ 

AM  =  Rl  ♦  RO  ' 

TMl  =  AM  *  (PTWDT(l)  -  DTWUTI 2 ) )  »  RHO  *  CP  »  DR _ 

PRINT  35,  AB,  AO,  AM,  TMI 

35  FORMAT I/6HAB  =  ,E12.4,6H  AO  =  ,FI2.4,5X,6H  AM  =  ,EI2.4,5X, 

I  14K  FIRST  TERM  =  ,612.4/) 

!  _ c _ 

HG  =  TMl  /  ( ( T0I2)  -  TOI II)  *  AB) 

!  TGAV  =  0.0  _ _  _ _ 

DO  50  I  =  1,2 

J  T_  =_  TOi  I ) _ _ _ 

TG  =  foil)  ♦  ( RHO*CP*DR*  AM  *  DTWDT 1 1 ) )  /  HG  /  AB 

1  _ TGAV  =  TGAV  *  TG _ 

PRINT  40,  HG,  TG 

40  FORMAT ( 5H  HG  =,E12.4,  IPX,  5H  TG  =,EI2.4) _ 

50  CONTINUE 

\  GAV  /  2. _ _ _ _ _ 

|  OTO  =  TO (2)  -  TOC  1 ) 

I _ TAVG  =  TO(l)  »  OTO  /  2. _ 

[  '  PRINT  55,  OTO,  TAVG,  TGAV 

:  55  FORMAT ( 18H  TEMP  DIFFERENCE  =,  E12.5,  15X,  23H  AVERAGE  OUTSIDE  TEMP 

T=,  E12.4,  IOX,  I4H  TGI  AVERAGE)  =,  E12.4/) 

_THE~rEfON'D~aUANTITY  IN  THE  FOLLOWING  STATEMENT  REPRESENTS  THE 
TOLERABLE  DIFFERENCE  BETWEEN  CONSECUTIVE  TEMPERATURES. 
FOLLOWING  NUMBER  (300.0)  MAY  HAVE  TO  BE  LOWERED 

IF<  DTP  .GT.  300.0)  GU  TO  100 _ _ _ 

I  C  REPRESENTS  THE  NUMBER  OF  ORDERED  PAIRS. 

_ 1 1  =  II  ♦  1 _ 

ID  =  II 

XIII)  =  TAVG _ 

Y  1 1  II)  =  HG 

Y2UI)  *  TGAV _ 

IM  =  II 

70  IM  *  IM  -  I _  _ 


IFIXI ID)  .GE.  X( IM)  )  GU  TO  80 
XD  =  X( ID) 


Y 2 ( 10)  »  Y2( IM) 
XI  IM)  »  XD 
Yll  IM)  *  Y 1 0 
Y2t  IMJ _«  Y2CL 
I D  »  ID  -I 
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oooo  o  oo 


g0  T0  70 

80  CONTINUE 

C  NOW  ORDERS 6  -  AND  II  VALUE  IS  THE  NUMBER  OF  ORDERED  PAIRS 

C 

100  CONTINUE 
C 

PRINT  105,  ixm,  nil),  Y2 (  I  )  »  I  =  1,  II) 

105  FORMATIIOX,  E12.5,  10X,  F12.5,  10X,  El?. 5) 

M  =  4 
NC  =  2 

CALL  LSTSOIM,  II, NC) 

110  CONTINUE 

111  CONTINUE 
CALL  EXIT 
END 


SUBROUTINE  LSTSQI M,N ,NC ) 

OIMENSION  B(20),C(420) ,D<420) 

COMMON  /BLK2/  X  <  30 ) ,  V1I3G),  Y2I10),  DTDTI30),  YYI30) 

THIS  SUBROUTINE  IS  GOOO  FOR  SECOND  OR  THIRD  DEGREE  POLYNOMIALS,  ONLY. 

M  =  NUMBER JjF  CONSTANTS  NEEDED  FOR  DEGREE  OF  POLYNOMIAL  DESIRED. _ 

"  . KEEP  THE  POL YN OM I AL  aT  DEGREE  TWO  TC  AVOID  POINTS  OF  INFLECTION. 

N  =  NUMBER  OF  ORDERED  PAIRS  READ  IN. 

DO  1  1=1, N 

1  YYI I )  =  Yll I ) 

NC  =  1  THE  FIRST  TIME  THROUGH  (MATCHES  TIME  VS  TEMP). 

NC  =  2  THE  SECOND,  AND  LAST,  TIME  THROUGH ( MATCHES  TIME  VS  HG  AND  TIME 
VS'TG(AVE)  ). 

DO  400  IJK  =  1 , NC 
N1=(M+1)*M 

DO  2  1  =  1, N1  _  _ _ 

2  D( I )=0. 

K =  1  _ _ _ 

3  XX  =  X ( K  ) 

Y  =  YY(K)  ‘ 

B  (  1 )  =  1 . 

DO  4  1=2, M 

I EXP= I-  1 

4  B(  I  )  =  XX  **  I E XP_  __  _  _ _ 

'1 1  =  0 

1  =  1 
J=1 
JJ=0 

5  I  SUB  1= 1 1  + 1 +J J 

I SUB2*  I  *J J  _  _ 

mSUBir=B(J)*B(ISUB2) 

0(  I  SUB  1  )  =  0(  I  SUB  1 )  *-C  (  I  SUB  1 ) 

IFU+JJ  .GE.  M )  GO  TU  6 
I-I  +  l 
GO  TO  5 

6  I«I+1 

fsuBi-imuj . 

ci tsuBi >«si j  )  ♦ y 

D( I  SUB  1 )  =  D( I  SUB  1 )  +C(  I  SUB  1 1 
1 F ( J  .GE.  M)  GO  TO  7 
II«II*MFl 

_ J  =  J  +  1  _  _ _ _  _ _ _ _ _  _ _  _ 

GO  TO  5 

7  K-K+l 
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1FIK  .LE.  N )  GU  TO  3 
J  J  =  M  + 1 
J=1 
1-2 

”k  =  0  .  '  " . . 

N  IN*  1  _ _ 

9  IJ0bl=JJ+J 

l J0B2* I +K  _ 

D< I JOB1 ) =D( IJOb2) 

IF(J_.GE. _NIN_)_GO  TO  11 _ 

. k-kvm*t'  ”  . 

j=j+i 

GO  TO  9 

11  N I N=N l N+  L  _ 

K  =  0 

I*!M  _ 

. j=i .  ' 

JJ f  JJ»M»1  _ _ 

IFljJ  .LE.  {M-l)*(M*m  GO  TO  9 

C  NORMAL  EQUATIONS  IN  Of  I )  COMPLETED  TO  HERE _ 

IN=  1 

1 OUT  =  M ♦ 1 _ _ 

DO  40  LP=1,M 

M8=M» 1 _  __  _ _ _ _ _ 

DO  20  I  =  1 , MB 

20  C(1)=D(1> _ _ 

IFIC(l)  .EC.  0.0)  GO  TO  800 

_ J  =  M»2 _ 

M  3= ( M-2 ) *M8 

K  =  0  _ _ _  _  _ _ 

25  DO  30  II*1,M 

I  IP* 1 1 *K  _  _ _ _ 

IIQ=J*II 

30  P(llP)  =  Q(IItt)-Q(J)*CIll*l)/C(l) _ 

d  1 1  iP*n*-o(  j)/c(  n 

J*  J_*M8  _ _ _ 

K=K+M8 

l  F  ( K  .LE.  M3)  GO  TO  25  __  __ _ 

DO  35  11*1, M 

_ I  IR*II»K _ 

35  D(UR)  =  Cm*l)/Cm 

DII!RHI>1./CI1I  _ _ _ 

40  CONTINUE 

JX  =  0  _  _ _ _ 

DO  50  !*l,Nl,M8 

_ JX*JX»1 _ 

B( JX)*D<  I  ) 

50  CONTINUE  _  _ 

PRINT  53 

53  FORMAT ( //53H  COEFFICIENTS  UF  FITTED  CURVE  BEGINNING  WITH  CONSTANT 
l//) 

_ PRINT  54,  ID! I ), 1*1, N1  ,M8) _ 

54  FORMAT  I  4020.8 ) 

IFIM  .  EU.  4  )  GO  TO  68  _  _ _ 

NBC* l *M0 

NAC»NUC*M8  _ _ 

8C«0(NBC) 

ACjlO  INAC) _ 

GO  TO  76 

68  NCC  «  l  ♦  M8  „  „ 

NBC  «NCC  ♦  MH 
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NAC  =  NBC  +  M8 

CC  =  D1NCC ) _ 

8C=D(NBC ) 

AC=P(NAC  > _ 

76  PRINT  77 

77  FORMAT I // 14X 1HX 19X1HY 1 7X8HF 1  TTEO  Y15X5HSLUPE ) _ 

IFIM  .NE .  3)  GO  TO  79 

_ BAA  =  -  BC  /  2.  /  AC _ ___ 

PRINT  78,  BAA 

78  FORMAT ( /25X,  29H  VERTEX  OF  PARABOLA  IS  AT  X  =  ,  E12.4/) _ 

GO  TO  84 

79  X1NFLP  =  -12.  »  BC)  /  3.  /  (3.  »  AC > 

PRINT  80,  XINFLP 

80  FORMAT  1 /IPX,  27H»»*  INFLECTION  POINT  AT  X  =,  E12.4,  4H  ***/> 

DISC  =  (2.  *BC )  **  2  -  4.  *  <3.  *  AC )  *  CC 

IFIDISC  »L  f»  0.0)  GO  TO  82 _ 

R1  =  (-2.  *  BC  ♦  SQRT(DISC))  /  2.  /  (3.  *  AC) 

R2  =  (-2.  »  BC  -  SQRT(OISC))  /  2.  /  (3.  ♦  AC) _ _ 

PRINT  81,  Kl,  R2 

81  FORMAT ( /28H  THE  REl  MAX  AND/OR  MIN  ARE  ,  5X,  F12.4,  IPX,  ^12.4/) 

GO  TO  84 

82  PRINT  83 _ 

83  FORMAT l /57H  THE  THIRD  DEGREE  POLY  HAS  NO  REL  MAX  NOR  REL  MIN  VALUE 

IS  /) _ 

IFINC  .NE.  2)  GO  TO  89 

84  IFUJK  .EQ.  2)  GO  TO  87 _ 

PRINT  86 

86  FORMAT I /15X,  34H  ♦  TEMP(OUTSIDE)  VS  HG  *»♦»♦/) _ 

GO  TO  89 

87  PRINT  88 

88  FORMAT  I / 15X,  34H  *****  TEMPI  OUT SI  DEI  VS  TG  *•***/) 

89  00  121  LLL-l.N _ 

XX  •  XtLLU 

Y  «  YYILLL  ) _ 

YF*0.0 

00  94  I»1,M _ 

IEXP-I-1 

IXI«XX _ _ _ 

IFIIXI  .NE.  0)  GO  TO  93 

_  IF! 1EXP.NE.  01  GO  TO  92 _ 

YF*YF*8( l ) 

GO  TO  94 

92  YF-YF 

GO  TO  94 

93  YF*YF*BI I)  *  XX  **1EXP 

_94  CONTINUE  _  _ _  _ 

IFIM  .EQ.  4)  GO  TO  96 
DT»2.*AC*XX  ♦  BC 
GO  TO  98 

C  96  OT  ■  3.  *  6(4)  «  XX**  2  ♦  2.  *  B(3)  *XX  ♦  8(2) 

96  DT  •  3.  *  AC  *  XX**'2  ♦  2.  *  BC  *XX  ♦  CC 

_98  IFINC  .EQ.  2)  GO  TO  99 _ _  _  _  _ 

IFINC  .EQ.  1)  YYILLL )  ■  YF 
DTDTILLL)  «  DT 
99  PRINT  120, XX,  Y,  YF,  OT 

120  FORMAT  I //4F20.6) 

121  CONTINUE 

IFINC  .EQ.  1)  GO  TO  800 _ _ 

IFUJK  .EQ.  2)  GO  TO  400 
DO  131  I  ■  1 »N 
lJl  YYI H  •  V2I I ) 
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400  CONTINUE 
800  RETURN 
END 


//PRINT  JOB  <0174093, II, 'CEO  STILES* 


JOB  16B 


C  ONE-DIMENS IONAL  TRANSIENT  HEAT  CONDUCTION  PROGRAM  IHT-2AI  00010 

C  THIS  PROGRAM  IS  A  GENERAL  PROGRAM  FOR  THE  SOLUTION  OF  CONDUCTION  00030 

C  PROBLEMS  WITH  TEN  OR  LESS  REGIUNS  INCLUDING  INTERFACIAL  RESISTANCES  00040 

C  BETWEEN  REGIONS  000*0  , 

OInENSION  lOFElOl 91  C 

OATA  I0F£t0  /*  • , • STEV* , *E  80* , • STWl • • *CX  • , • AMSW* , *E-RE • , • T-E  C 

|A*,*  •/  C 

WRITE i 14)  IOFELO  C 

DIMENSION  XLAR(*),VLA8<*>,GLAB<S>,DATLAR<S)  C 

DIMENSION  XBST  S I, VBSI SI »GRS|  SI ,OBS<  5)  C 

DIMENSION  LABL  1  (  201  ,LASL2<  201  ,L  ABL  3420)»LABL4I20)  .LABLSUOI  .LABLAI  C 

120) • L ABL 7<20)«L AREA) 20)  C 

DIMENSION  LABI 9 1 20) ,LA6L 101 20) ,LA0ll II 20) .LABL 12 <201 • LABL I  3120 1, l A  C 

IBL  14  <  20 )  ,L  ABL  IS<  20)  ,LABL  16120)  «LAOl  17120) , LA<<LIR<20I  , LABL  19120) , LA  C 

2 BL 20 1  20) • LABL 2 1120)  C 

DIMENSION  LABL 221 20) ,LABL23< 20) »L ABL 24 1 20) , LABL 2S < 201 , LABL 26 1 20) «l  C 

IABL27<  20 ), LABL 2B( 201  * LABL 291 20) *LABL 101 20) , LABL 3 1 120)  C 

READ  S00,XLAB,YLA8«GLAB,0ATLA8  C 

READ  SOO.LABL l,LABL2,LA8L3,LABL4,LABLS,LABL6,LA6L7,LABLB  C 

REAO  SOO, l ABL 9, LABL 10,1 ABL I 1 ,L ABU 2 »L ABL 13, l ABL 14, LABL IS, LABL I 6, LA  C 

IBL 1 7, LABL 18, LABL 19, LABL 20, LABL 2 1  C 

REAO  S00,XSS,YBS,GBS,0BS  C 

REAO  S00,LABL22,IA8L23,LABL24,LARL2S,LABL26,LABL27,LABL2R,LABL29,L  C 

1  ABL  30,1  ABL  31  C 

SOO  FORMAT <20A4 I  C 

01  ME  NS  ION  TRA0I20I,  RA0ISTI20) 

DIMENSION  TORI  100) « 70S 1 1001 • TNEI 100)  C 

C  ENO  PLOT  OIMENSION  C 

DIMENSION  TTItSO),  RH0PI20I, 

1  MCI 20 1 •  TG<20I*  TTII20) ,  CPII20I 

C  00070 

^•DEFINITION  OF  LABELEO  COMMON  ~  BLRI.BLR2,  ANO  SIRS  00080 

COMMON  /BLR1/  T< ISOI »C < ISO) ,CXI ISO! ,H| ISO) »MXl ISOI , IBOOVI 10,2 >  00090 

COMMON  /BLR2/  RAOI II l I ),NOO(S< 101 , XX/120) ,BE TAI 10) »CP<20) »RH0I20I* 
2fMIS$,RM02,CP<, XXR2»B0VR< III, RID SOI, Rill  ISO) ,DR< |0),AI9),|TB<I|)  001 10 

C  001 10 

C*4|N|T IALI2ATIUN  OF  VARIABLES  NOT  10CATE0  IN  LABELED  COMMON  00140 

OATA  RRUN,  TNUN. TDENOM.Of ,OT|NEX,DOTX, I X.NBOOV/  I* 

2  >0,  t*«  l«,  »00S,  . 2S,  3,  2/ 

c  00170 

C**REAO  CHARACTERISTICS  OF  PROBLEM  —  RAW  INPUT  OATA  00180 

C  00190 

CR40EFINIT ION  OF  NAN  ANO  NAN|  00200 

NAMELIST  /NAM/  T,RRUN, T NUN, TOE NON,  WOOES*  XRl,BETA,CP»AHO.BOtA, 

2EM|SS*0<*  0t|MEX,0DTX,U,NB00Y.C*2,RH02,XRR2.RA0«l,A,ITB,NR0,NSC 

4/NAM | /OT |MEX,00TX,0f ,11 , NBOOY, I X, XXR/,RHOI,CPl ,EM| SS* TNMN, (Of NUN,  00240 

S  NOOES*  A*|!B*WAO,MSC 

C  N2N»2  CHECKS  NO.  OF  ROUNOS  AS  A  MULTIPLE  OF  2. 
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H2H»2 
KKfta  0  i 

C*a*a*aa**aNCPLUT  «  0*  -I.  OR  1  •••••••••••••••*••♦••••••••••••••••••••••••••••• 

(•••••••••a  NCPlOT  «  -I  WHEN  NO  PIO!  IS  OESIRfO  . . . 

(•••••••••a  nCPLOT  •  0  WHEN  BOTH  KAO|Al  AMD  ACS  PLOTS  ARE  OF SIRE O*************** 

C#*a*aa***a  NCPLOT  ■  1  WHEN  ONLY  UNF  Of  THf  TnQ  TYPES  OF  PLOTS  IS  OESIREO******* 

Cim*m**mn|RPU  a  o»  OR  1  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
C, »••••*••••  NBRPLT  a  o  WHFN  THE  RAOIAL  OR  BOTH  THE  PLOTS  ARE  Of  SIRED* ••••••••••• 

C*aaa******  NBRPLT  »  I  WHEN  ONLY  THE  ACS  PLOT  IS  OESIRfO************************ 

C**********WHEN  NCPlOT  a  o,  NBRPLT  HOST  ALSO  a  0.  ••••••••••••••••••♦•••••••••*- 

NCPLOI  «- l 
NBRPLT  •  0 

If  I  NCPLOT  .EQ.  -I  »  CO  TO  2 
If  I  NBRPLT  .  EO  .  I  I  CO  TO  *60 
NRAO  *  0 

IF  I  NCPLOT  .EO.  0  )  Ct)  TO  *60  < 

CO  TO  2 
*60  K8RSK  >  | 

TBRIK6RSKI  a  TO. 0 

TOSIKARSK)  >  70.0  i 

TMEIKBRSK)  a  0.0 
2  RCXOIS.NAMI 

c 

C  FOLLOWING  CAROS  ARE  REOUIRFO  IF  OESIRE  TO  RETAIN  T|WF  PER  BURST • 

C  BUT  A  LOCATION  CHANGF  REQUIRES  VARYING  BARREL  THICKNESS. 

C  020LD  IS  ORIGINAL  02  VALUF  IN  FEET. 

020L0  •  .02TSB 

C  02NEH  IS  PRESENT  VALUE  OF  02. 

OINEN  *  02 

OriNEA  a  020L0  ••  2.  *  DTlwf*  /  P2NFW  ••  ? 

,C 

C 

IOUT  *  0 

00  )  I  •  laNBOOY 
1  I OUT  a  IOUT  •  NOOfSill 
IUUT  a  IOUT  *1 

TOUTS  •  moufi 

CALL  TGHGI TOUTS. TOTS* AVGhG I 

TUI  a  TOTS 

•OVRttl  •  HR 2  /  AVGHC 
C 

C  IF  NCR  STAYS  IFRO  THEN  THERE  IS  AN  EVEN  NO.  0*  ROUNDS  PER  BURST. 

NCR*0 

|«<*  SO/NRO 
IFIIM.CQ.  01  I RR* I 
IFIN00INR0.N2NI  .10.  01  GO  TO  * 

NCR*  t 

OtlNfRa  OIINE 1/2. 

c  IRR  REPRESENTS  THf  BURST  TO  Af  PRINTED  (EVERY  IAX  BuRStl. 

C  NAD  REPRESENTS  THf  NO.  OF  ROUNDS  PER  BURST. 

C  NSC  REPRESENTS  THE  NO.  OF  SECONDS  OF  COOLING. 

*  PRINT  2D2.NRO.NSC 

202  FURWAIIlHl.ISR.IS.IBM  AOUNO  BURSTS  AMO  SECONDS  COOLING/ I 

Nf»»NA0/2 

C  NIB  IS  NUMBER  OF  PASSES  PER  BURST 

IPINCK  .EO.  II  NIB  «  NRO 
C  NlC*NlB*NSC*10/ I  (IF  A  MULTIPLE  OF  »l 
NIC  •  NIB  •  NSC  •  10  /  2 

C  NIC  IS  NUMBER  uf  passes  per  BURST  PLUS  NUMBER  Thru  coining 

■  NNIC  •  INK  -  NIRI  /  2  *  NIB 
NANaO 
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o  o 


00270 

00280 


c 

C**CALCULATE  DIMENSIONLESS  LUMPED  PARAMETERS.  HX(I)  AND  C ( f  > 
C 

CALL  LUMP  ( I I.NBODV* OZ. TT1.CP1 I 


C**WRITE  PROBLEM  PARAMETERS 
WRITE(6,NAM1) 

WRITE! 6. 5 1 

5  FOP "AT I 7H0REG10N. 3X5HIB0DY* 3X  9HRADI I CFTI .5X6HDRCFTI ,5X8HB0YR(FT) , 

26X2HCP. 8X3HRHO. 8X2HKZ , 6X4HBETA  } 

WRITE (6. 7)  I J. I BODY! J. 1 ) , (BODY! J « 2 ) .RADI I(J).OR(JJ ,BDYR(J) .CPI J) • 
2RH0! J ) . XKZ I J I .  BETA(J),  J*l, NBODY) 

7  FORMAT (I3,I8,I4,3E12.3«  F 10. 3, 2F 10. 1 ,F 11.6) 

I  *  NBODY  ♦  1 

WRITE (6*  9) I. RADI  1(1) .BDYR ( I ) 

9  FORMAT! 13, 12X.E 12. 3. 12X.E12.3//J 

WRITEI6, 11) 

11  FORMAT (3XIHI, 7X5H  H( I),  12X4HCU),  12X4HT(I),  7X6HRA0IUS  ) 
WRITEI6. 13) ( 1.  H( I ) , C < I ),T(1),RI(! )«  1  =  1,11) 

13  FORMAT (I4,2E16.4,F13.2,F13.5) 

JI  =  1 1- 1 
DO  90  1=2, Jl 
RACIST ( 1-1)  =  Kill) 

90  CONTINUE 
JJJ  =  1-1 
C 

C**CALCULATE  OR  INITIALIZE  VARIOUS  QUANTITIES  -  SAVE  T(I)  AND  OTIMEX 

TSEC  =  0Z**2*RH0Z»CPZ*3600./XKRZ 
IIM1  =11-1 
I IM2  =  II  -  2 
IIPl  =11+1 
DO  15  1=1, IIPl 
15  rT ( 1 1  =  T!I) 

ATIME  =  DTIHEX 

DDDTX=DDTX 

N=0 

NMN  =  0 
TAUT  =  .0 

CALL  TAVE! II, I IP1) 

C 

SQIN  =  0.0 
SQSTR  =  0.0 
SQUT  =  0.0 
NT  =  0 

DO  17  J  *  1, NBODY 
NT  l  =  NT  ♦  l 
NT  =  NODES! J)  ♦  NT 
00  17  I  *  NTl.NT 

17  RHOP(l)  *  RHO! J  ) 

C 

C 

C**$T ART  OF  SOLUTION  OF  PROBLEM 
C 

18  CONTINUE 
TOUTS  =  T(IOUT) 

CALL  TGHG! TOUTS, TOTS, AVGHG) 

T ( 1 )  =  TOTS 

BDYR! 1 )  »  XKRZ  /  AVGHG 

CALL  LUMP  ( 1 1, NBODY, DZ, TT1,CP1 ) 

C 


00300 

00310 

00360 

00370 

003BO 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 

00490 


00500 

00510 

00520 

00530 

00540 

00550 

00610 

00620 

00630 

00640 

00660 

00780 


00820 

OOB30 


LEECH 
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00840 


C  POINT  OF  MAJOR  LOOP  ENTRY  —  SN25(N0  NEW  OTIMEX),  SN24INEW  OTIMEX) 

IFURET  .EQ.  2)  00  TO  25 
00  19  1*2, I  INI 

19  CXN)  *  Cm/0T1NEX*2.  00860 

25  NNMM  *  NMN 

CALL  CHANGE  I NBODY, T SEC, TAUT *11*1 X,N,NMN,NIB»NIC ) 

NMN  s  NNMM 

CALL  SOLVE  (IIMl, IIN2. II, NBODY, BETA, TAUT) 

C 

DETM  *  OTIHEX  *  TSEC  /  3600.0 
QSTR  =  0.0 
00  27  I  «  1*  NT 

QSTR  »  (Rill  +  I)  -  Rlim*  RHQP(I)  *  CPU)  *  6.2832  *  RIM)  *  IT 

III)  -  TT I  I))  /  DETM  ♦  QSTR 
27  CONTINUE 

QIN  HXm*XKRZ*6.2832*{T{l>  -  TI2) ) 

QOUT  x  HXI I INI )*XKRZ*6« 2832*1 T< I INI )  -  Till)) 

SQIN  =  SQIN  ♦  QIN  *  DETM 
SQSTR  *  SQSTR  *  QSTR  *  DETM 
SQUT  =  SQUT  ♦  QOUT  *  OETM 
SUN  *  SQSTR  *  SQUT 
ENBL  =  (SQIN  -  SUN)  /  SQIN  *  100.0 
C 

DO  29  I»1,IIP1 
29  TTII)  *  Til) 


N  *  N  ♦  1  00890 

NMN  ~  NMN  ♦  1 

TAUT  =  TAUT  ♦OTIMEX  00900 

C  00910 

C**IF  N/IX  IS  AN  INTEGER  CALCULATE  VARIOUS  QUANTITIES  AND  STORE  ITENS  00920 

CALL  TAVElIItUPl)  00960 

C  01010 

C**END  OF  TIME  STEP  01020 


IF (NMN  .NE.  NIB)  GO  TO  31 
OTIMEX  *  OTIMEX  ♦  2. 

WRITEI6*  203)  OTIMEX, NIB 

203  FORMAT!  /21H  NMN=NIB  AND  OTIMEX  *,E12.5,10X,6H  NIB  =  ,I5> 

GO  TO  33 

31  IFINHN  .NE.  NIC)  GO  TO  33 
OTIMEX  =  OTIMEX  /  2. 

WRITE (6, 204)  OTIMEX, NIC 

204  FORMAT (  /21H  NMN-NIC  AND  OTIMEX  -,E12.5,10X,6H  NIC  =, 1 5 ) 

33  IFINCK  .NE.  1)  GO  TO  35 

IRET  »  1 
GO  TO  37 
35  IRET  =  2 

37  IF  I NMN  .EQ.  NIC)  GO  TO  41 

IFIMOOINMN,  25)  .NE.  0)  GO  TO  18 
GO  TO  42 

41  NMN  =  0 
NXM=NXM*l 

42  IFIMQDINXM, IXX )  .NE.  0)  GO  TO  18 

CALL  RESULT(TAUT,IIM1,II,TNUM,T0EN0M,DZ,NB0DY, QSTR, DETM, SQSTR, ENBL 
2, IIP1, SQIN, SQUT, TIME) 

IF  (  NCPLOT  .EQ.  -I  )  GO  TO  331 

IF  (  NCPLOT  .EQ.  1  .AND.  NBRPLT  .EQ.  1  )  GO  TO  321 

KRAD  =  KRAD+1 

TME(KRAD)  *  TIME 

DO  311  1*2,15 

TRADU-l)  *  Til) 

311  CONTINUE 
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CAL).  GRAPH)  jJJ,  RACIST,  TRAD,  11,1,  10.5,  10. 5,0. 0 .O.O.XL AB, YLAB.GL AB, 0 
IATLAAI 

ALL  LE UFA  1 2.0, 9.8, 0.2.LABL1 .0.80.0,0,0.0.0.0,0,0.0,0) 

CALL  LETTER)  1.0, 1.6. 0.1. LA8L 2, 0,80,0,0,0,0,0.0.0 .0*0.0) 

CALL  LETTER)  1.0, 1. A, 0.1.LA8L8, 0,80, 0,0, 0,0,0, 0,0, 0,0,0) 

CALL  LETTER) I. 0, 1.2,0. l.LABL A, 0,80, 0.0, 0,0, 0,0,0, 0,0,0) 

CALL  LETTER) 1.0, 1.0,0. l.LABL 3,0, 80.0, 0,0,0, 0.0,0. 0.0,0) 

CALL  LET  TER) 1.0, 0. 8, 0.1.LABL6, 0,80, 0,0, 0,0,0, 0,0, 0,0,0) 

GU  TO  1700, 710. 720,710. 740,750, 760,170), KRAD 

700  CALL  LETT6R11. 0,0. 6.0. I, LABL5, 0.80, 0,0,0, 0,0,0, 0,0, 0,01 
GU  10  BBS 

710  CALL  LETTER)  1.0,0. 6,0.1, LABL9, 0.80. 0.0. 0.0. 0.0. 0.0, 0.0) 

GO  TO  880 

720  CALL  LETTER) 1.0,0. 6,0.1, LABL1C.0, 80, 0,0, 0,0, 0,0, 0,0, 0.0) 

GO  TO  888 

730  CALL  LETTER)  1.0, 0.6,0. l.LABLU, 0.80. 0,0, 0.0, 0,0, 0,0, 0.0) 

GO  TO  888 

740  CALL  LETTER) 1.0. 0.6,0. 1.LABL12, 0,80,0, 0.0, 0,0, 0,0, 0,0.0) 

GO  TO  888 

760  CALL  LETTER) 1.0. 0.6, 0.1, LABLl 3. 0.80. 0.0,0. 0,0, 0.0. 0,0. 01 
GO  TO  888 

760  CALL  LETTER) 1.0. 0.6.0. 1.LABL14, 0,80, 0,0, 0,0, 0.0. 0,0, 0,0) 

GO  TO  888 

770  CALL  LETTER) 1.0. 0.6, 0.1, LA8L15, 0,80,0, 0.0. 0.0. 0.0, 0,0. 01 

888  CALL  LETTERll. 0, 0.2, O.l.LABL 7, 0,80, 0,0, 0,0,0, 0,0, 0,0,0) 

IF  1  NCPLOT  .EQ.  0  I  GO  TO  321 
GO  TO  331 

321  KBRSK  •  KBRSKM 
TBRtKBRSK)  -  T ( 2 ) 

TOS(KBRSK)  «  T I IUUT) 

TMEIKBRSK)  •  TIME 

C  NXM  INDICATES  WHICH  BURST  HAS  BEEN  FIRED. 

331  IFINXM  .GE.  1)  GO  TO  45 

C  ABOVF  IF  STATEMENT  INDICATES  NUMBER  OF  BURSTS  TO  BE  FIRED. 

GU  TO  18 

C  01200 

C**RESET  INITIAL  CONDITION  AND  TIME  INCREMENT  —  READ  NEXT  CASE  —  SN26  01210 
45  KKR  ■  KKR  ♦  1 

IF  !  NCPLOT  .EQ.  -I  I  GO  TO  365 

IF  (  NCPLOT  .EQ  .  0  .AND.  NBRPLT  .EQ.  0  )  GO  TO  360 
IF  (  NBRPLT  .EQ.  0  )  GO  TO  365 

360  CALL  GRAPH(KBRSK,TME,TBR,11,1,10.5,10. 5 ,0,0, 0,0 ,XBS, YRS.GBS .DBS) 

CALL  GR APH) KBRSK, TME.TOS, 11, 1,0,0,0,0,0,0,0,0,0,01 
CALL  LE T TER ( 2. 0, 9.5,0. 2,LABL22, 0.80. 0.0.0, 0.0, 0.0, 0,0,0) 

CALL  LETT ER 12.0,9. 1,0.2,LA8L23,0,80, 0,0, 0,0, 0,0,0, 0,0,0) 

CALL  LETTER(2.5,2.3,0.1,LABL24,0,80,0,0,0,0,0,0,0,0,0.0) 

CALL  LE T TER ( 2.5. 2. L.O.l.LABL 25,0, 80,0, 0,0, 0,0, 0,0, 0,0,0) 

CALL  LET  TER  I  2. 5, 1.9,0. 1, LABL26, 0.80. 0.0. 0.0, 0,0. 0,0. 0,0) 

CALL  LET  TER  1 2. 5, 1 . 7, 0. 1 ,LA8L2 7,0, 80,0, 0,0, 0,0, 0,0,0, 0 ,0 » 

CALL  LETTER) 2.5, 1.5,0. l.LABL 28,0, 80.0,0, 0.0, 0,0, 0,0, 0,0) 

CALL  LET  TER (2. 5, 1. 3, 0. 1 .LA8L29.0, 80,0,0, 0,0, 0, 0,0, 0,0,0 ) 

CALL  LETTERI2.5, 1.1,0. l.LABL 30,0, 80,0, 0,0, 0,0, 0,0, 0,0,0) 

CALL  LETTER(2.5,0.9,0.1,LABL31,0,80,0,0,0,0,0,0,0,0,0,0) 

365  IFIKRUN  .GE.  KKR)  GO  TO  370 
DTIMEX  «  ATIME 
00  47  1 -  1, I  IP  1 
47  Tl I )  ■  TT)  I  > 

GO  TO  2 

370  CALL  EXIT 

END  01260 

SUBROUTINE  LUMP) 1 1 ,NBOO Y, 02 , TT1 ,CP 1 ) 
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dimension  ttusoi,  rhdpi/o). 


1  HGI20),  TGI20),  TT1I10),  CP1I10) 

c 

common  /eu  1/  T  (  ISO), Cl  ISO)  .CXI160I.HI  |  SO  t  ,MX(  ISO)  ,  IB0DYI10.2I  012B0 

COMMON  /BLK?/  ft  *1)1 1 1  1 1  )  .NODE  SI  l  U )  .  XK  Z  I  20 »  ,BE  TA I  1 0 )  .C  P I  20 )  ,RmCI20)  , 

2EM 1 SS, RHQZ  ,CPZ»XKRZ»Bl)YR(  11 1  ,ft  1  (  ISO  I  .ft  I  1  II  SO)  ,0«  (  1 0»  01300 

C  01310 

C**T  M I S  SUBROUTINE  CALCULATES  THE  DIMENSIONLESS  LUMPED  PARAMETERS  01320 

AZ=RHOZ*CP2*0Z**2  01330 

Cl  1  .0  013*0 

cm  =  .o  013S0 

I F I HDYR {  1 ) . EO. . 0  1  GO  TO  3  01360 

HKI1I  >  RADI  111 l/BOYRI 1  )  01370 

Ml  1 1  »  MX1 1 1  01 JHO 

I  BODY  I  1,1)  «  2  01390 

GO  TO  S  01*00 

3  I  BODY  1 1 , l)  »  1  01*10 

5  Rill)  »  RADI  1(1)  01*20 

C  01*30 

C**BEGlNNING  OF  LOOP  TO  CALCULATE  C1I)  AND  HID  FOR  NROUY  REGIONS! J)  01**0 

DO  9  J* l • NBODY  01*SO 

ORR  ■  RAC  1 1 1 J»  1 )  -  RADI  II J )  01*60 

DKIJ)  >  ORR /FLOAT  I  NODE  SI J ) -  1 )  01*70 

I  BODY ( J , 2 )  *  I  BUOY  I J , 1 )  ♦  NODE  S ( J )  -  1  01*80 

IB  *  I  BODY  I J, I )  01*90 

IE  «  180CVI J , 2 )  -  1  01600 

All  IB)  «  RA01I ( J)  01610 

C  01620 

C**CALCULAT ION  OF  Cl  I  I  AND  Hi  I  I  FOR  REGION  J  01630 

00  1  I  *  I B, ! E  01670 


C 

C**USE  FOLLOWING  IF  WANT  VARIABLE  CP. 

C  1 X 1  •  1 

C  TEMPI  *  TUI 

C  CALL  LINEARIT6MP1,TT1,CPI,CPJ, I  XI ) 

C  CPU)  a  CPJ 

C 

IF  I J  .NE.  2)  CPU)  *  .11 
I F I J  .EO.  2)  CPI  1  I  -  .IB 
AJ  *  ft HO  I J ) • CP  I  I )*0R I J I  /AZ 

CUB)  «  AJ*  IRK  IB  )  ♦  DR  I  J  )  /*.  )  /2  .  ♦  C!  01660 

CALL  XKKSII.XKZJ) 

XK 2  I  I  )  «  XK  Z  J 
I F I J  .EO.  2)  XKZI  1 1-1.0 
XKDR  -  XKRZ  •  DR  I J ) 

BJ  «  XK Z 1  I ) /(XKRZ PORI J) ) 

C 


HU)  »  BJ* IRl(l)+DR(J)/2.)  016B0 

R 1 1 1  *  1 )  ■  R  I II  )  *  DR  I J  )  01690 

1  CUM)  «  AJ*R  1 1  l  ♦  1 1  01600 

Cl  I E ♦ 1 )  •  AJ*(RIIIE*l)-DR(J)/*.)/2.  01610 

C  01620 

C**CECK  TO  SEE  IF  INTERFACIAL  RESISTANCE  IS  ZERO  AND  PROCEED  ACCORDINGLY  01630 
IFIBDYRI J*1 I.EO..O)  GO  TO  2  016*0 

Cl  ■  .0  01660 

I  BODY  I J* !• 1 )  *  (BODY! J, 2 )  *  1  01660 

HXIIEM  )■  R 1 1  I E ♦  1  ) /BDYR I  J*  1 )  01670 

HI  IE  *1)  a  HXI 1 E m 1 )  01680 

GO  TO  9  01690 

2  Cl  ■  Cl IE*1)  01700 

I  BODY  I J* 1, 1 )  >  I  BODY  I J, 2 )  01710 
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9  CONTINUE  01/20 

IFIBOVRI NBOUY* l ) . NE . . 0 )  GO  TO  tl  01/10 

II  ■  If  ♦  l  01/4" 

go  rn  l i  oi / so 

11  II  -  Ifc  *2  01/60 

cm  i  =  .o  ni  //o 

13  M( 111  =  .0  01 /HO 

RIUI)*  RAO  I  1 1 NBODV  ♦  1)  01/00 

C  01  MOO 

C**C ALCUL  A I  £  t  HE  DIMENSIONLESS  RAOIUS  Rll  01810 

00  16  IMill  01820 

16  RIMI)  -  (Rid)  -  RADII) 1))/1R AD  lit  NBOD  V*  1 )  -  RA01 1(1)1  01810 

RETURN  01840 

END  OlflSO 

SUBROUTINE  SOLVE  ( 1 1  Ml , 1 1  M2 , I  1 , NROO V , BE T A , T AU T I 

DIMENSION  GE( 1 50  I  *  F  E ( 1  SO  I »  DE ( 150)»BETAI 101 ,BE  I ISO) •  BII1S0)  018/0 

COMMON  /BLK 1/  T ( ISO) »C( ISO) ,CX( ISO) *H( ISO) «HXI ISO) .1  BODY ( 10,?  )  OlHHo 

C  01890 

C**CORRECT  THE  BOOV  CONDUCTANCES  FOR  VARIABLE  CONDUCTIVITIES  01/00 

1  DO  3  J*  l *NBODV  01910 

IB  =  I  BODY ( J, 1)  0192" 

I  E*  I  BODY  (  J  ,  2  t  -  1  01  MO 

DO  3  1*18, IE  01940 

3  HX ( I  1  =  HU»*fl.  ♦  BETA(J)*)T(I)  ♦  TIIU))/2.)  019S0 

C  01960 

C**START  OF  ELIMINATION  --  C (ANK-NICOL  SON  ALGORITHM  019/0 

00  9  1  =  2,  I  IM1  u 1 9 rt o 

Cl  =  HX(  II  ♦  MX)  1- 1  I  0  1990 

BF  (  I  1  =  cxm  ♦  Cl  02000 

9  ei m  *  cxm  -  ci  02010 

GE I  2 )  *  BE ( 2 )  02020 

F  E ) 2 )  *  (BI)2)*T(2)  ♦  H X ( 2  » • T «  3 »  ♦  HX(1)*T(1)*2.I/GFI2)  020  10 

DO  5  I*3,IIMl  02040 

DEIII  *  -HX) I-1)/GE< 1-1 1  02050 

GE1I)  *  BE  <  I  )  ♦  HX<|-n*DE(l)  02060 

5  FEII1  =  IHXIIIMIIMI  ♦  MX(  l-l)»Tf  1-1)  ♦  BI<I>*TI|>  ♦  HXU-ll*  020/r. 

2  FE ( |-1I)/GEI I  I  U20H0 

FE(IIMl)  «  FEIIIMll  ♦  HX) I  I  Ml  )  *T( II l/GEIIIMl)  02090 

C  02100 

C**BACK  SUBSTITUTION  02110 

T(IIMl)  *  FEIIIMll  02120 

DO  7  1*2,  1  IM2  02.1  10 

J  *  II  -  I  02140 

7  T  ( J I  «  F  E  ( J 1  -  OE(J*ll*T( J*l)  021S0 

RETURN  02160 

END  021/0 

SUBROUTINE  RESULTITAUT.IIMl, I  I , TNUM, TDENCM, 02 , NBOOV, OSTR , DE TM , 

2  SOSTR.ENBL,  I 1P1, SQIN, SOUT, T TME I 

DIMENSION  TSTAR) 1501, XM< 101,  V(SOO) 

COMMON  /BLK 1/  T ( ISO) «C( ISO) , C  X ( ISO) i H I  ISO) ,HX( ISO) , IBODY ( 10,2)  02210 

COMMON  /BLK2/  RADIK  111, NOOES) 10) , XK2 t 20 ) ,BE TA I 101 ,CP) 20) ,RHO( 20  I , 
2EMISS,RH02,CP2,XKR2,B0YR( lll.RK 1 50 ) • R I  I  1 150) , DR (  10)  02230 

C  02250 

C**CALCULATE  DIMENSIONAL  TIME ,HFAT  FLOWS  PER  UNIT  DEPTH,  TSTARS,  M'S  02260 

C**Af,0  WEIGHTED  AVERAGE  TEMPERATURE.  PRINT  THESE  QUANTITIES.  022/0 

CALL  TAVEI 1 1 • I  IP  1 )  022H0 

TSEC  «  02**2  *  RHD2*  CPZ  •  3600./  XKR2  02290 

TIME  ■  TAUT  •  TSEC  02300 

Q1N  *  HX( 1)*XKR2*6. 2832*1 Tl 1)  -  T(2))  02310 

QOUT  *  HXI  IIMl)*XKR2*6.2832*milMl)  -  TUI))  02)20 

PRINT  100,  SQIN,  SQSTR,  SQUT,  ENBL 

0 
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100  FORMAT ( 5X,  |)H  OJN  SUM  *  ,E12.*»  5X,  15H  QSTORED  SUM  »  .  E12.*»5X, 
t  12H  OOUT  SUM  »  ,  E 12. A/  3SX,  20H  ENERGY  BALANCE  t  •  .  EI2.*) 
PRINT  31, TSEC. TAUT, TIME 

31  FORMAT  I / 7H  TSEC  -.E12.5.9H  TAUT  *,F12.5.9H  TIME  «,E12.5/I 
00  1  I-1.IIP1 

I  TSTARII)  »  ITU!  -  TNUMI/TOENOM 
WRI TE (6, 5 )  TAUT 

5  FORMAT  I ////22HO  DIMENSIONLESS  TIME  ■  , F  7. 3  , 10X28HHt AT  FLOW  PER  FT  ( 

2BTU/HR-FTII 

WRITEI6, 7)  TIME, OIN, OOUT 

7  FORMAT ( 22H  REAL  TIME  I  SECONDS  I  * ,E 1 1. 3, 3X*HQI N*E 1 2. 3, 7h  Q0UT*.E12 

2.31 
C 

C**PR INT  THE  DIMENSIONAL  TEMPERATURES 

WRITE (6, 111  Tl  II, IIMl.l T( II, I >2, II  Ml | 

II  FORMAT  I  /35H0THE  DIMENSIONAL  TEMPERATURES  ARE  /6H  TI1)=,F10.2/ 
213H  T 121  THRU  T I » I  3, 9H I  FOLLOW/ I 5f 10. 2 . 5X, 5F 1 0. 2 )) 

WRITEI6, 13)11,  Till),  TIIIP1I 
13  FORMAT ( 3H  T I , I  3, 2HI * ,F 12. 2, 6X, 7HT I AVE I » , F 1 2. 2 ) 

RETURN 

END 

SUBROUTINE  TAVEI I  I , I  IP  1 ) 

COMMON  /BLK 1/  T( ISO  I, Cl  ISO l, CXI  150 ) ,H 1 1 SO ) , HX 1 1  SO  I , I  BODY  I  10,2) 


C 

C**CALCULATE  WE IGHTFD  AVERAGED  TEMPERATURE  AND  STORF  IT  IN  TIIIP1I 
SUM  >  .0 
SUM2  »  .0 
00  39  |>l,ll 
SUM  •  SUM  ♦  Cl  I  I • T C I  1 
39  SUM2  >  SUM2  ♦  cm 
T I  I  IP  1 )  a  SUM /SUM 2 
RETURN 
END 

SUBROUTINE  CHANGE  INBOOY, TSEC , TAUT , 1 1 , 1X,NNN,NMN,NI R , NI C ) 

DIMENSION  HZ  I  111, Nit 11)*N2( 11) 

COMMON  /BLK 1/  T I  1501 ,C ( 150) ,CX I  150 » .HI  150) ,HX < 150) , I  BODY ( 10 , 2 ) 
COMMON  /8LK2/  RADI  II 1 1 1 ,NOOES< 10) , XK» 1 20) ,BE T A 1 10 ) .CPI  20 ) , RHC < 20 ) , 
2EMISS,RH0Z,CPZ,XKRZ.B0YRI 111 ,R1( 1501 ,RI 1  1 150) .OR  I  10) , A<9) , I TBI 11 ) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


J  a  NUMBER  OF  R'S  WHICH  ARE  TEMP.  OR  TIME  DEPENDENT 

Nl(J)a  RESISTOR  NUMBER  --  Nil J)  ■  Jl 

N2CJI-  RESISTOR  TYPE 

HZ  I J )■  RESISTOR'S  INITIAL  VALUE 

A  >  ARRAY  CONTAINING  COEFFICIENTS  FOR  FUNCTIONS.  EXPONENTS  ETC. 
TSEC  «  CONVERSION  FACTOR  (REAL  TIME  IN  SECONDS  »  T|MEX#TStCI 
EXPOla  EXPONENT  N  WHERE  H  ■  HZ'ABSITIJl)  -  T I Jl ♦ 1 ) I **EXP01 
|TB  >  ARRAY  CONTAINING  TYPE  KEY  FOR  ALL  BOUNDARY  RESISTORS 
TYPE  -I  Ha  CONSTANT 

TYPE  a  2  Ha  HZ*F3f TIME) 

TYPE  «  3  H«  HZ9(0T|9*EXP01 

TYPE  ■  A  H  a  hr  ♦  HC 

TYPE  a  5  H  a  HZ»F5(T|ME)  —  F5  IS  A  PERODIC  RECTANGULAR  WAVE 

STORE  INITIAL  VALUES  AND  DETERMINE  WHICH  RESISTORS  ARE  NOT  OF  TYPE  1 
Ntei«NIB«l 
N1CI-NICM 

IFINMN.EO.  NIC)  NMN»0 
IFITAUT.GT..O)  GO  TO  1 
N-0 


023*0 
02370 
02  360 
02  390 
02*00 
02*10 
02*20 
02*50 
02*60 
02*70 
02*60 
02*90 
02500 
02510 

02680 

02690 

02900 

02910 

02920 

02930 

029*0 

02950 

02960 

02970 

02980 

02990 

03000 

03020 

03030 

03050 

03060 

03070 

03080 

03090 

03100 

03110 

03120 

03130 

031*0 

03150 

03160 

03170 

03180 

03190 

03200 

03210 


03220 

03230 


NMNaO 

N8  •  IF  I XI  At •)  )  032*0 
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S4  =  imi  4i4)  i 
T 1  -  1(1) 

Til  =  till) 

im  --  |  -  1  SS*  .  I  714t  -R/xyr/  a. 

E  XPU1  «  A  (  7 1  r0f 

J  =  0 

lHIIrMM.EY.il  L.!  I  10  I 
J  =  1 

SHI)  ‘  l 
s/t  i  i  =  nmn 
h  /  (  n  *  h  x  1 1 1 

7  C( .  s  1=1,  Shi. i  r 

IMIlnilMl.iu.il  O’)  lo  •> 

J  =  J  *  1 

SlIJl  =  l  HOt*  V  I  1  ,  /  ) 

S  /  l  I  )  =  I  I  H  l  1  »  L  I 
J  l  :  .  1  (  J  I 

0/ ( J 1  =  HXtjll 

S  CoSIlSUl 

c 

C • • PO I  NT  UE  ENTRY  EUR  T I MF  .G 1 .  7  E Rl!  --  CAlClllAlf  Nr  W  t»  0  V  Tt  MPfRATURFS 

C 

C 

1  (IMF  =  T  Al)T  •  T  SF  C 

Till)  ■  T  I  1  •  (  l .  ♦  TIMf*(A(5l  ♦  A(4)*T|ME)) 

C 

C*  *  I  F  ALL  R*S  ARE  COS  S I  AS  I S  RETURN  OTHERWISE  RECALCULATE  THOSE  CHANGING 
IMj.lu.1l  RE  TURN 
U(J  11  IM.J 
J  l  =  (.HI) 

OTEMR  *  A8S  (  TIJM-Tljni)) 
ir  (GIEMP.EU..O)  OTEMPsl. 

M  *  N  /  (  1  I 

r.n  to  ( 1 1,  i/.  l ),  ia,  lsi  ,m 

1/  HX(JI)  =  H/lll  *11.  ♦  A(S)*SIN(  A(M*T|Hf  |  » 

GO  TO  11 

l)  IUIJ1I  *  HE  (  |  )  •  UTfcMP  *♦!  *P()1 

GO  TO  U 

14  TA  *  T  (  J  1  I  ♦  4,60. 

T  8  =  I  ( J 1 ♦  11  ♦  S60 . 

HXIJ1)  =  HR  *  R I  I J 1  I  *(TA**/  ♦  TB*  > • ( T A  ♦  IB » 

/  ♦  H/lll  •  CUSP  EXPUl 

GO  TO  l l 

is  ieinmn.lt.  NlMt  Hxijin  0/(1) 

IMNMN.GT.  N|R  .ANO.  NVN.LT.  Nl  C 1  I  HX(J1)  »  M/m  •  A  (  S  ) 

IF  (N.IM.N'TI  S  =-| 

S  «  N  ♦  l 
NMS  =  NMN»  I 
11  CONTINUE 

SSI  =  NSN  ♦  l 

IE ( (MOO (NS  l, IX 1 .NF.O) .OR. I J.EQ.O)  1  RE  TURN 
PU  ?\  1 « 1 . J 
Jl  *  Nl  I  l I 

21  T(iniM)  >  HXIJI)  •  XKR2  /  RI(Jl) 

HfcTuRN 

ESO 

SUBROUTINE  XKK S ( I , XX  1 

COMMON  /flLRl/  T< ISO) ,CI ISO) ,CXI 1S0I ,M( ISO) ,HX( ISO) , 180071 10.2) 

T  T  « T  (  I  ) 

1MTT-IA72.)  10,10,14 
10  XA-/8. )0-.00870*TT 


0  5/SP 
0  5/60 

U  5/80 
0  5/40 
0)500 
05)10 
0  5  5/0 
0  )  )10 
0  5  54,0 
0  5  5S0 
0  5  560 
0  5  570 
0  5  580 
05  540 
0  5**00 
05410 
054/0 
0  54  50 
05440 
0  54SO 


05480 
0  5440 
O5S0O 
05S  l  O 
05S/0 
05S50 
0  5S40 
05SS0 
0  5S60 

05SI50 
0)*>90 
0  5600 
0)610 
056/0 
056  50 
03640 
056SO 
0)660 


0)690 

0)700 

0)740 
037SO 
0)760 
0)770 
05780 
0)790 
0  )A00 
03810 
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14 

20 

JO 


1 

c 

10 

c 

100 

2 

1 


c,  U  111  2.) 

XK  =  10.  1')*.0014  7*11 

CUNT  1 NUl 

RETURN 

END 

SUBROUTINE  L INEARI A.X, Y.VV, l ) 

DIMENSION  X I  20 ) «  Y I  20 ) 

IFIYI  !♦!  »  .LT.  Ylll  I  f.O  TO  too 

USE  EOLLUM ING  IF  AS  V  INCREASES  X  INCREASES 

ifi  A-am  i  j.2,2 

USE  FOLLOWING  IF  AS  Y  INCREASE  S  X  OECRFASfS 
IFIA-XC I  II  2,2, J 
1  =  1*1 


*0,. 


GO  TO  l 
1=1-1 

vv=vin*(A-xii«in/<xin-xii«in«Y(i»n*iA-x(ni/ixii*n-xiin 


RETURN 
f  Nil 

SuflROuT INE  TGHGl TOUTS. TOTS.AVGHG) 

ITS  =  TOUTS 

C  RATIO  =  1.268 

IFITTS  .LT.  1S9.2)  TTS  =  1S9.2 

AVGHG  *  • 8669S6E- 0 J  •  TTS  ••  2  -  .622TT  •  TTS  ♦  312.44 

TOTS  =  .1161224F-01  •  TTS  ••  2  -  .4496S  •  TTS  ♦  1606. IS 

AVGHG  >  AVGHG  *  .20  *  AVGHG 

TOTS  *  TUTS  ♦  .40  •  TOTS 

RETURN 

ENU 

BLUCK  OATA 
C 

C**INIT IAL I  Z AT ION  OF  LABELED  COMMON  TO  DEFAULT  VALUES 

COMMON  /BLR  1/  T( ISO) .Cl 1S0I .CXI ISO).H(ISO) .HXI ISO) , IBODYI 10.2) 


C 

CUMMON  /BLR  2/  RA01M 11). NODE  SI  10) . XR2 I  20) .BE T A  I  1 0 ) ,CP I  20 ) .RMU I  20 ) . 


C 

2EMISS.RH02.CRZ.XKRZ, BOYR I  1 1  I .R 1 1  1  SO ) ,R i I < 1  SO > ,DR I  1 0) . A| 9 ) , I T B I  1 1 ) 
DATA  EM|  SS.RHOZ.CRZ . XRRZ.XRZ.BE T A, CP. RMU/ 

C 

2  1.,  490. ,.11,  10..  20*10.,  10*. 0,20*. 11.20*490.0 

C 

)  /.NODES. T, BOYR/ 10*S.  1 . , 149* . 0, l l *. 0/ , A, I TB/6*. 0, . 2S .2* .0 , 1 1 • 1 / 
END 


01B20 
010  JO 
03H40 
01BS0 


0  TRIO 


03910 

01920 
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